Bose-Einstein condensate as a quantum memory for a photonic polarization qubit 
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A scheme based on electromagnetically induced transparency is used to store light in a Bose- 
Einstein condensate. In this process, a photonic polarization qubit is stored in atomic Zeeman 
states. The performance of the storage process is characterized and optimized. The average process 
fidelity is 1.000 ± 0.004. For long storage times, temporal fluctuations of the magnetic field reduce 
this value, yielding a lifetime of the fidelity of (1.1 ± 0.2) ms. The write-read efficiency of the pulse 
energy can reach 0.53 ± 0.05. 
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I. INTRODUCTION 

Optical quantum-memories [l| based on electromag- 
netically induced transparency (EIT) Q are a very active 
research area. Specifically, storage of classical light pulses 
and of single photons 6,, 7} was demonstrated. In 
addition, different polarizations of light were stored ei- 
ther using extended atomic level schemes or by 
converting the polarization into other degrees of free- 
dom before storage [Iol - fl3l |. The vast majority of appli- 
cations envisioned for quantum memories requires that 
quantum entanglement is first generated between two or 
more particles and that the quantum states of one or 
several of these particles are subsequently stored in a 
quantum memory. The crucial point is that this entan- 
glement must survive the storage. Recently, this aspect 
was experimentally demonstrated in three independent 
experiments [14| - [l6j . 

Here we report in detail on the performance of the 
quantum memory used in one of these experiments [l6j . 
The experiment uses an 87 Rb Bose-Einstein condensate 
(BEC) to realize a quantum memory for the polarization 
qubit of a single photon. A Raman scheme based on EIT 
is used to implement storage and retrieval of the photon. 
The atomic level-scheme is extended to allow for stor- 
age of the photonic polarization qubit in two atomic spin 
states. Quantum process tomography is used to deter- 
mine the process fidelity which quantifies how well the 
polarization is maintained during storage. In addition, 
the decay of the process fidelity with increasing storage 
time is monitored. 

All experiments reported here use classical light pulses 
instead of single photons, thus profiting from count rates 
which arc much higher than in Ref . [1 6j j . These increased 
count rates yiel d a more precise value of the process fi- 
delity. Ref. [16[ reports that storage and retrieval cause 
no discernible deterioration of the fidelity of the entan- 
gled state within an error bar of several percent. The 
present experiment still observes a process fidelity com- 



patible with unity, but now with an error bar that is 
an order of magnitude smaller, thus demonstrating more 
clearly the capabilities of the BEC as a quantum memory. 

In Ref. [lfjj ]. a single 87 Rb atom in an optical high- 
finesse cavity was used to generate a triggered single pho- 
ton in such a way that the photon's polarization qubit is 
entangled with the spin state of the single atom. This un- 
paralleled source is combined with the BEC that serves 
as a quantum memory for the single photon. The BEC is 
well suited for this purpose because the absence of ther- 
mal motion allows for long storage times, the large optical 
depth allows for high write-read efficiencies, and excellent 
internal-state preparation allows for high-fidelity storage 
of a qubit in atomic spin states. The experiment thus 
combines two different systems, each ideally suited for 
its purpose. The resulting hybrid character of the sys- 
tem poses an experimental challenge because the dipole 
traps that hold the single atom and the BEC in place have 
depths of several millikelvin and several microkelvin, re- 
spectively. Due to the resulting ac-Stark shifts, the single 
photons generated from the single atom are blue detuned 
by typically 70 MHz relative to the free-space atomic res- 
onance, whereas the ac-Stark shifts experienced by the 
BEC are negligible on this scale. In the present paper, 
we therefore experimentally study the efficiency of light 
storage in the regime of 70 MHz single-photon detun- 
ing. The regime of much larger single-photon detunings 
has been studied theoretically |17H20| and experimentally 
[l3|,[2l[ before, but those results are not immediately ap- 
plicable to our system. 

The paper is structured as follows: Section HI1 describes 
the experimental implementation, Sec. lIIII shows how well 
the polarization is maintained during storage, and Sec. 
II VI studies the write-read efficiency. Appendix [Al presents 
a simple model for coarsely estimating the write-read ef- 
ficiency. 



II. EXPERIMENTAL IMPLEMENTATION 



A. Electromagnetically induced transparency 
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EIT employs a control light beam to manipulate the 
propagation of a probe light beam inside a medium. 
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FIG. 1: Atomic level scheme of the Di line in 87 Rb. Probe 
light (dashed arrows) with an arbitrary superposition of the 
polarizations a + and a~ couples the initial population (•) 
in the hyperfine ground state \F, tti_f) = ] 1, 0) to the Di-line 
excited states \F',m' F ) = |1,±1). 7r-polarized control light 
(solid arrows) couples these states to the hyperfine ground 
states \F, thf) = |2, ±1). This makes it possible to store the 
probe-light polarization qubit in the qubit space spanned by 
the atomic states \F,mF) = |2, ±1). 



Light storage in EIT-based schemes relies on the fact that 
the group velocity of the probe light t> gr can be reduced 
compared to the vacuum speed of light c by many orders 
of magnitude [13] by choosing a small value for the con- 
trol intensity. Upon entering the medium, the temporal 
duration of the probe pulse remains unchanged, whereas 
its spatial length is drastically reduced due to the small 
group velocity. A pulse which in vacuum is much longer 
than the medium can thus be fully compressed into the 
medium. 

Once the pulse is fully inside the medium, one can 
ramp the control intensity to zero in an adiabatic fash- 
ion. In our experiment, we implement an approximately 
linear, temporal ramp of the control intensity which lasts 
30 ns. This is sufficiently adiabatic according to Ref. [23| . 
This ramp reduces v gI all the way to zero and the pulse 
is stopped inside the medium. The pulse is stored for 
a time i s torc which can be chosen freely. After this, we 
ramp the control intensity back on and the pulse resumes 
its propagation [3H51 l24j. During the storage time t storc , 
the pulse exists in the medium in the form of an atomic 
spin wave. If the light pulse is compressed such that 
it fits inside the medium, then the spin wave stores the 
longitudinal and transverse shape of the light pulse. 



B. Atomic level scheme 

Figure [T] shows the atomic level scheme used in our 
experiment. Control and probe light for EIT are both 
resonant with the atomic D-y line of 87 Rb at a wavelength 
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FIG. 2: Simplified scheme of the optical beam path, 
detailed description is given in the text. 
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of A = 795 nm. The atoms are initially prepared in the 
hyperfine ground state \F, mp ) = 1 1, 0) . The a ± polar- 
ized components of the probe light couple this population 
to the excited hyperfine states \F',m' F ) = |1,±1). The 
7r-polarized control light transfers this population to the 
ground hyperfine states \F,mp) — |2,±1). 



C. Optical beam path 

Figure [5] shows a simplified scheme of the optical beam 
path. An 87 Rb BEC serves as a quantum memory. The 
BEC is illuminated by 7r-polarized control light propagat- 
ing along the y axis with a waist (1/e 2 radius of intensity) 
of ~ 100 /im which is much larger than the Thomas-Fermi 
radii of the BEC so that the control light intensity can be 
approximated as constant across the BEC. In addition, 
the BEC is illuminated by probe light propagating along 
the z axis focused to a waist of 8 /im. This is comparable 
to the Thomas-Fermi radii so that the probe beam sam- 
ples some fraction of the transverse inhomogeneity of the 
BEC. To obtain a well-defined transverse mode for the 
probe light before impinging onto the BEC, the light is 
sent through a single-mode fiber. The polarization of the 
probe light can be er + , er~, or any superposition thereof. 

After storing and retrieving the probe light, we need 
to measure its polarization. To this end, the beam path 
ends with a polarizing beam splitter (PBS) cube and two 
identical detectors, one in each output port of the PBS. A 
quarter-wave plate (QWP) followed by a half-wave plate 
(HWP), both placed right in front of the PBS allow for 
the selection of an arbitrary polarization basis. As detec- 
tors, we use photomultiplier tubes (PMTs) in this paper, 
instead of the avalanche photodiodes that we used in Ref. 
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D. Stray light filtering 

Stray light is an issue in our setup. Much of it is elimi- 
nated using mechanical shielding and temporal gating of 
the detector signals. The remaining stray light is dom- 
inated by control light off-resonantly scattered from the 
BEC during the retrieval of the probe pulse. 

This stray light level would be unproblematic for the 
measurements presented here, but the experiments re- 
ported in Ref. 16] required a substantial suppression. 
In the beam path from the BEC to the detectors, our 
setup therefore includes a single-mode optical fiber for 
transverse mode filtering and a filter cavity for spectral 
filtering. 

The spatial filtering with the single-mode fiber makes 
use of the fact that storage and retrieval have little effect 
on the transverse mode of the probe light, whereas the 
control light is off-resonantly scattered from all positions 
in the BEC and into all directions. The fiber reduces the 
stray light power that reaches the detector by a factor 
of 0.068. In the absence of the BEC, the fiber reduces 
the probe light power by a factor of 0.66. Hence, the 
fiber increases the signal-to-background ratio by one or- 
der of magnitude. Storage and retrieval in the BEC com- 
promise the transverse mode of the probe beam slightly. 
This causes an additional reduction of the probe light 
power by a factor of 0.88 at the single-photon resonance 
and by a factor of 0.80 for a single-photon detuning of 
A c = 2ir x 70 MHz. In addition, the fiber suppresses the 
excitation of higher transverse modes in the subsequent 
filter cavity. 

The filter cavity is a near-planar, plane-concave Fabry- 
Perot resonator with a finesse of 180 and a free spectral 
range of 40 GHz. The cavity has a transmission of 0.8 
at the resonance, which is tuned to the probe light fre- 
quency We expect that the scattered EIT control light is 
either elastically scattered or Raman scattered, thereby 
transferring an atom from F = 1 to F = 2. Hence, we 
expect the scattered light to be red detuned from the 
probe light by either 6.8 GHz or 13.6 GHz. Transmission 
through the cavity suppresses the light power for each 
of these frequencies by a factor of 2 x 10 -4 . The cavity 
length is stabilized against long-term drift with a piezo 
actuator using a Pound-Drever-Hall technique with light 
at a wavelength of 830 nm. This light is overlapped with 
the probe light on a dichroic mirror in front of the mode- 
filtering fiber. The 830 nm light transmitted through the 
filter cavity is separated from the probe light using an- 
other dichroic mirror behind the filter cavity. A small 
fraction of the 830 nm light keeps propagating towards 
the detectors. This light is removed with dielectric inter- 
ference filters which are not shown in Fig. [2] In addition, 
the 830 nm light source is turned off during any time 
intervals where relevant detector signals are expected. 

The optical fiber for stray-light filtering is not polar- 



ization maintaining because it must work equally well 
for all possible input polarizations. Hence, the polar- 
ization after transmission through the fiber is related to 
the input polarization by a unitary transformation. As 
long as the fiber is not moved mechanically or exposed to 
temperature changes, this transformation is temporally 
stable. The same applies to the polarization transforma- 
tions caused by the filter cavity, by the dichroic mirrors, 
and by other mirrors in the beam path which are not 
shown in Fig. [5J We compensate the resulting overall 
transformation using two QWPs followed by a HWP - 
a combination which can generate any unitary transfor- 
mation. 



E. Preparation of the BEC 

We produce an almost pure BEC in the hyperfine state 
\F,mF) = |1, — 1), using radio-frequency (rf) induced 
evaporative cooling in a magnetic trap, as described in 
Ref. [25[ . The gas is transferred into a crossed-beam op- 
tical dipole trap operated at a wavelength of 1064 nm. 
The measured trap frequencies are (ui x , uj v , uj z )/2it = 
(70, 20, 20) Hz with gravity pointing along the x axis. 
A magnetic hold field of ~ 1 G applied along the z axis 
preserves the spin orientation of the atoms. 

We use two consecutive microwave pulses, each with 
a pulse area of tt, to transfer the population into the 
internal state needed for our EIT level scheme. Starting 
from state |1, — 1), the first pulse transfers the population 
into state |2, 0). Subsequently, the second pulse transfers 
the population to state |1,0). The total process trans- 
fers ~ 90 % of the atoms into state |1,0). Atoms left 
in the F = 2 hyperfine states are then removed with 
blast light. This is followed by temporary application 
of a strong magnetic field gradient which removes atoms 
with mp from the shallow optical dipole trap. Af- 
ter this procedure, the total atom number in undesired 
internal states lies below the detection limit of our setup 
which we estimate to be ~ 200 atoms. For the rest of the 
experiment, the magnetic hold field is reduced to typi- 
cally 0.1 G. At this point, the BEC typically contains 
N = 1.2 x 10 6 atoms. The corresponding Thomas- Fermi 
radii are (R x ,R y ,R z ) = (7,25,25) fj,m. 



III. A QUANTUM MEMORY FOR THE 
POLARIZATION QUBIT 

We now study how well the polarization of the probe 
light is maintained during storage. The probe beam prop- 
agates along the z axis. Hence, an arbitrary incoming 
polarization state can be expanded as c+|<r + ) + c_|<r~) 
with coefficients c+ and c_ . With the atomic level scheme 
shown in Fig. [T] this state is mapped onto the atomic 
state c+|2, +1) + c_|2, —1). The retrieval process maps 
the atomic state back to the original polarization state. 
In our experiment, both mapping processes work ex- 
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FIG. 3: Faraday rotation. The normalized Stokes parameter Si /So oscillates as a function of storage time due to an applied 
magnetic hold field. The line is a fit of Eq. ([5]). The best-fit value for the e -1 ^ 2 damping time is a a = (1.1 ± 0.2) ms. Note the 
breaks on the horizontal axis. 



tremely well, but magnetic field noise causes a deteri- 
oration of the state for long storage times. 

We apply a magnetic hold field of B z ~ 0.1 G and 
orient it along the z axis. This suppresses undesired 
transitions between different Zeeman states caused by 
components of the magnetic-field noise perpendicular to 
the z axis. In Sec. lIII Al we discuss the Faraday rotation 
caused by this hold field. In Sec. IIII Bl we use quantum 
state tomography to characterize the deterioration of the 
polarization for long storage times and show how tech- 
niques that reduce the magnetic-field noise improve the 
performance of the system. 



A. Faraday rotation 

The polarization of the probe light can be character- 
ized using the Stokes parameters [2g, l27| 



So — Ih + Iv > 
S2 = Id — I a, 



Si — Ih — Iv, 
S3 = Ir — II , 



(la) 
(lb) 



where Ih, Iv, ■ ■ ■ denote the intensity detected after a po- 
larizer that transmits only one polarization, namely hor- 
izontal H, vertical V, diagonal (+45°) D, anti-diagonal 
(—45°) A, right circular R, or left circular L. Here, R 
and L correspond to a + and er~ . The Stokes parameters 
can be regarded as the components of a four-dimensional 
Stokes vector. So describes the total intensity, whereas 
the three-dimensional vector 



u 




(2) 



describes the polarization of the light and is well suited 
for graphical visualization, in close analogy to the Bloch 
vector. We call u the Poincare vector. Its unit sphere is 
called Poincare sphere. 



The magnetic hold field along the z axis gives rise to a 
Faraday rotation of the Poincare vector around its z axis 
at an angular frequency 



uj f = B z 



(3) 



where [is = 27rfi, x 1.40 MHz/G is the Bohr magneton 
and g F is the Lande factor. For the levels used in our 
experiment gp = 1/2 and Amy = 2. The total rotation 
angle of the Poincare vector 



(j) = W_F(t sto rc + Td{L)) 



(4) 



has one contribution t s t ore from the storage time with the 
control light off and another contribution Td{L) which is 
the delay of the probe pulse caused by the propagation 
through the medium of length L with the control light 
on. Note that the rotation of the polarization vector of 
the electric field in real space is a factor of 2 slower than 
the rotation of the Poincare vector. 

Figure [3] shows experimental data of this Faraday ro- 
tation for a linear input polarization. The line shows a 
fit of a sinusoid with Gaussian damping 



Si 
So 



e -*?tare/2<^ COS(0~ <po), 



(5) 



where a a is the e" 1 / 2 damping time, where <j> is given by 
Eq. (jj), and where <po represents the input polarization. 
The best-fit values are ujf = 27r x 0.20 MHz, yielding 
B z = 0.14 G, and a a = (1.1 ± 0.2) ms. 

A careful inspection of the experimental data points for 
long times in Fig.[3]leads us to an interesting observation, 
revealing the physical origin of the damping. Unlike the 
best-fit curve, the data points, which represent a single 
experimental shot each, do not show a decrease of the 
peak-to-peak amplitude. Instead, they are noisy insofar 
as they do not form a smooth sinusoid. We attribute 
this to irreproducible, temporal fluctuations of B z , which 
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effectively scatter the data points along the horizontal 
axis. Only if we were to average several experimental 
shots to represent their mean value, would we observe a 
reduction of the pcak-to-peak values in the experimental 
data. 

A fit of Si /So = Acoa(4>-4>o) to the data in Fig.[3]with 
^storc < 20 /is yields a best-fit value of A = 1.02 ± 0.04. 
The fact that A is consistent with 1 shows that tempo- 
ral fluctuations of B z on a time scale of 20 /is or faster 
have no discernible effect. Temporal fluctuations of B z on 
slower time scales manifest themselves in our experiment 
only as shot-to-shot fluctuations of <fr — (/i^c^Ami? /K) 
J* st °" !+Td ^ L '> dtB z (t). In our experiment, each shot re- 
quires the preparation of a new BEC, which takes 20 s. 
This suggests that shot-to-shot fluctuations of B z prob- 
ably yield an important contribution to the shot-to-shot 
fluctuations of <ft. 

B. Quantum process tomography 
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FIG. 4: (Color online) The damping factor a of Eq. © as 
a function of storage time. The lowest damping is obtained 
when synchronizing the start of the EIT write-read cycle with 
the 50 Hz ac line voltage (■). Alternatively, we can use an 
open-loop feed-forward compensation (•). Without any com- 
pensation, the damping is much stronger (a). The lines show 
Gaussian fits according to Eq. (J5J). 



We now turn to a complete characterization of the ef- 
fect which the storage and retrieval process has on the 
polarization. As the process does not conserve the to- 
tal intensity, a full description of the process must use 
the four-dimensional Stokes vector, not just the three- 
dimensional Poincare vector. We consider the regime 
of small probe intensity. Hence, the dependence of the 
Stokes parameters S° nt of the retrieved probe pulse on 
the Stokes parameters Sf 1 of the incoming probe pulse 
can be approximated as linear 

X ( 6 ) 

3 

M is called Miiller matrix [H,[27j]. 

As shown in Fig. [21 we use a PBS with detectors behind 
both output ports. According to Eq. ([1]), a full character- 
ization of the Stokes vector requires such measurements 
for 3 different settings of the wave plates in front of the 
PBS, which select the measurement basis. This set of 
3 measurements fully characterizes the quantum state of 
the polarization and it can be regarded as quantum state 
tomography [28]. 

To determine M experimentally, we use a set of 4 lin- 
early independent input Stokes vectors (e.g., H, D, R, 
and L) and perform quantum state tomography of the 
output state generated for each input state. This set of 
12 measurements allows for a complete determination of 
M and it can be regarded as quantum process tomogra- 
phyH. 

Performing such quantum process tomography, we find 
that in our experiment the Miiller matrix is always well 
approximated by 

/ 1 \ 

yr _ acoscj) — asin0 j ,„x 
^ I asincf) a cos ' ^ ' 

\ 1 / 



where </> is the angle resulting from the Faraday rotation, 
a is a damping factor, and 77 is the write-rcad efficiency 
which is experimentally found to be independent of the 
input polarization. 

Above, we concluded from Fig. [3] that the polariza- 
tion at short storage times is essentially pure, whereas 
at long storage times shot-to-shot fluctuation of </> must 
be taken into account. Hence, each individual shot can 
be described by some realization of M as in Eq. ([7]) 
with a = 1 and with some value of (j) which exhibits 
shot-to-shot fluctuations. We assume that the values of 
(j) have a Gaussian distribution with root-mean-square 
(rms) width u^. Averaging over many shots yields Eq. 
© with a = exp(-<r|/2). 

To develop a simple model for the dependence of 
a on tstorc, we assume that only shot-to-shot fluctu- 
ations of B z contribute to the fluctuations of </>, i.e. 
we approximate B z as constant during each individual 
shot. With this approximation, the values of B z will 
have a Gaussian distribution with rms width <tb and 
with (70 — (psgF Aitif /fo)o~Bt store, where we neglected 
T d (L) < i s torc- This yields 

a = cxp(-i s 2 torc /2 ( r2) (8) 

with 

1 duj F 
a a 6B Z 

For a linearly polarized input state, Eqs. (J7]) and <j8j) re- 
produce Eq. ([5|) . The Faraday rotation is a unitary evolu- 
tion. For any given i s torc, its effect can be compensated, 
e.g., using wave plates and it is therefore not much of 
a concern. The non-unitary damping a, however, irre- 
versibly deteriorates the performance of the memory. 

Experimental results for the time dependence of the 
damping parameter a, as determined by quantum pro- 
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cess tomography, are shown in Fig.|4j Data taken without 
any reduction of magnetic- field noise (A) yield a value of 
a a = 0.06 ms for the e -1 / 2 damping time, corresponding 
to c b = 2 mG according to Eq. @. In our experiment, 
the majority of this magnetic-field noise is periodic and 
in phase with the 50 Hz ac line voltage. We can sup- 
press this noise drastically by synchronizing the start of 
the EIT write-read cycle with the ac line voltage. The 
corresponding data (■) in Fig. [4] yield a best-fit value 
of a a — (1.0 ± 0.1) ms, corresponding to &b = 0.1 mG. 
Evidently, the synchronization improves a a by a factor 
of ~ 20. The data in Fig. [3] were also taken with this 
synchronization and essentially reproduce the improved 
value of a a . 

Our experiments described in Ref. [l6j required a rep- 
etition of EIT write-read cycles at a rate of 10 kHz for a 
total time span of several seconds. Hence, write-read cy- 
cles had to occur at essentially all possible phases of the 
50 Hz ac line voltage. To reduce the noise in these mea- 
surements, we first determined the values of B z {t) for one 
50 Hz period in a series of calibration measurements. We 
then ran a current through a coil to compensate the re- 
corded noise with an open-loop feed-forward circuit. The 
corresponding data (•) in Fig. 2] yield a best-fit value of 
<j a = (0.49 ± 0.04) ms, corresponding to <jb = 0.2 mG. 
This compensation was good enough not to be the lim- 
iting factor in the overall experiment of Ref. , where 
we observed the same lifetime but with an error bar that 
was a factor of four larger. 

The full information from the quantum process to- 
mography is contained in M. To compare the overall 
performance of different quantum memories, one often 
uses the average process fidelity as a figure of merit. In 
terms of quantum states, the fidelity can be written as 
F = Tr(pi n/ o out ), where Tr denotes the trace, p denotes 
the density matrix, and we assumed that p m represents a 
pure state. For polarization states, this can be rewritten 
as F = (1 + u in ■ ii out )/2. When averaging this quantity 
over all possible pure input states, i.e. over the surface 
of the Poincare sphere, we obtain the average process fi- 
delity (F) 0, El. 



rotation, i.e. for 



After compensation of the Faraday 
: 0, Eq. © yields 
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FIG. 5: Write-read efficiency r\ vs. Rabi frequency of the 
control light fi c . All data were taken at the two-photon res- 
onance. Data taken on the single-photon resonance A c = 
(o) display a clear maximum as a function of Q, c - Data taken 
at a single-photon detuning of A c = 2ir x 70 MHz (•) display 
a lower maximum efficiency. The lines are fits of the simple 
model from appendix A to the data. 



read cycle 



Er, 



V 



Eh, 



(11) 



defined as the energy of the retrieved probe pulse -E re tr di- 
vided by the energy of the incoming probe pulse E- m . For 
our memory, rj is independent of the polarization of the 
probe field, as seen in Eq. ([7]). Hence, for understanding 
rj it suffices to consider the case where the probe polar- 
ization is fixed to a + . This simplifies the relevant atomic 
level scheme to a A-type three-level system. Unlike pre- 
vious experiments by other groups, our work has a focus 
on the regime of 70 MHz detuning from the single-photon 
resonance. 

We experimentally study the dependence of rj on the 
intensity of the control laser in Sec. lIV Al In Sec. lIVBl we 
study the decay of the efficiency for long storage times. 



(F) = i(2 + a). (10) 

The synchronized data in Fig. [4] yield (F) = 1.000 ± 
0.004 at t stoIC = 1 us and (F) = 0.90 ± 0.02 at i st0 re = 
800 /is. The value at t B tore = 1 M s shows that the state 
mapping between photonic and atomic qubit states works 
extremely well. 



IV. WRITE-READ EFFICIENCY 

Now, we turn to another important figure of merit for 
light storage, namely the efficiency of a complete write- 



A. Dependence on the control intensity 

Figure O shows the experimentally observed depen- 
dence of rj on the Rabi frequency of the control light 
f2 c . Data were taken for a storage time of 1 fj,s with 
an incoming probe pulse that has a Gaussian intensity 
profile 

I iD {t,z) = I exp(-^(t-^yy (12) 

Here, Jo is the peak intensity and r p is the temporal rms 
width of the intensity. Data in Fig. [5] were taken for 
T p = 94 ns. The control light was turned off to = 230 ns 
after the maximum probe intensity entered the medium. 
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Data taken at the single-photon resonance A c = (o) 
show a maximum of 77 at O c ~ 2ir x 20 MHz. This value 
of Q c agrees fairly well with the prediction of the sim- 
ple model developed in appendix A. The observed max- 
imum efficiency of 77 ~ 30 %, however, is a factor of ~ 2 
lower than the expectation from the simple model. We 
attribute this discrepancy to the simplicity of the model 
and to experimental issues, such as inaccuracies in the de- 
termination of the experimental parameters. Note that 
because of the results of Fig. [5l the data in Figs. [3] and 0] 
were taken at f2 c = 2ir x 20 MHz, where 77 is maximized. 

We find experimentally that the write-read efficiency 
77 is increased if we slowly decrease the intensity of the 
control beam while the probe pulse enters the medium. 
This observation agrees with a more rigorous optimiza- 
tion of 77 for a homogeneous medium in Ref. [301 ] . Hence, 
all data in Fig. [5] were taken with a linear ramp of the 
intensity of the control beam applied, with a ramp speed 
that is experimentally found to maximize 7;. This ramp 
is not included in our simple model. For storage, the 
control light is on for 300 ns because prior to this, there 
is no probe light inside the medium. The horizontal axis 
in the figure shows the Rabi frequency corresponding to 
the time-averaged value of the control intensity during 
this 300 ns control pulse. 

In a measurement independent from Fig. [5j we 
achieved a write-read efficiency of 77 = (53 ± 5) %. The 
gain in efficiency compared to the data in Fig. [5] resulted 
from two changes in the experimental procedure. First, 
we removed the filter cavity and the mode-filtering fiber 
after the BEC and, second, we truncated the Gaussian 
input probe pulse in time, such that it misses exactly that 
part of its falling edge that cannot be stored anyway be- 
cause it reaches the BEC after f2 c is already ramped to 
zero. 

For reasons discussed in the introduction, our experi- 
ments in Ref. [16| had to be operated at a single-photon 
detuning of A c = 2ir x 70 MHz. An investigation of 77 
at this detuning was therefore necessary. Experimental 
results (•) are shown in Fig.[5j These data display a max- 
imum value of 77 - 20 % at fi c ~ 2ir x 20 MHz. The value 
of Q c at which the maximum occurs is essentially identi- 
cal to the data at the single-photon resonance, whereas 
the maximum efficiency is further reduced. The physical 
origin of this reduction is discussed in appendix IA 31 

The lines in Fig. [5] show fits to the data, based on the 
simple model developed in appendix A. More precisely, 
the dash-dotted line from Fig. [5] was taken and two fit 
parameters were introduced, each representing a linear 
scaling, one for 77 and one for f2 c . 



B. Dependence on the storage time 

Our experiments in Ref. [l6| also required an investiga- 
tion of the time scale on which 77 decays during storage. 
Thermal motion is known to be the limiting physical ef- 
fect in many experiments. Using a BEC or an optical 



lattice, however, thermal motion can be suppressed dras- 
tically, resulting in a very slow decay of 77(i s tore) [3~TI — 
33]. Unlike those experiments, our experiment does not 
use co-propagating probe and control beams. Instead, 
the level scheme shown in Fig. [T] requires the two beams 
to propagate perpendicularly to each other. The result- 
ing differential photon recoil is much larger than for co- 
propagating beams. In our experiment, the lifetime of 
7?(^store) is predominantly limited by this recoil, similar 
to Ref. Q. 

As a result of the photon recoil, atoms in hypcrfme 
states F = 1 and F — 2 move relative to each other. If 
after t store these two atomic clouds do not overlap any 
more, the retrieval does not produce a directed beam 
and hardly any signal reaches the detector. In our exper- 
iment, the control and probe beams propagate along the 
y and z axes, respectively. Hence, the differential pho- 
ton recoil incurred in the Raman transition is directed in 
the yz plane, where the BEC is symmetric with Thomas- 
Fermi radii R y — R z = 25 [atxi. 

The single-mode fiber between the BEC and the de- 
tector poses an additional constraint, also related to the 
photon recoil. Not only do the two atomic clouds need 
to overlap, the emitted light must also match the trans- 
verse mode of the single-mode fiber, resulting in a spatial 
filtering in the xy plane. The fiber is pretty well mode 
matched to the incoming probe beam, which has a beam 
waist of w = 8 /im, thus setting a length scale for the 
spatial filtering that is more stringent along y than the 
Thomas- Fermi radius R y . 

To obtain a simple estimate for 77(i s tore)j we approxi- 
mate the medium as homogeneous, which is justified by 
w <C R y and w R z . In this approximation, only the 
recoil along the y axis is relevant, i.e. the photon recoil 
of the probe laser is irrelevant. As the mode of the opti- 
cal fiber has a Gaussian transverse profile, the decay of 
??(istorc) is expected to be Gaussian 

i?(W«>) = 77(0)e-*--/2^ (13) 

with a e -1 / 2 time = raw / \/2hk c = 1.0 ms, where m is 
the atomic mass and k c = uj c /c is the wave vector of the 
control beam. Note that the factor v2 here has nothing 
to do with the modulus of the differential photon recoil. 

Figure [5] shows experimental data for 77(i st oro), re- 
corded for an essentially pure BEC and a single-photon 
detuning of A c = 2tt x 70 MHz. The line shows a fit of 
Eq. (fl3|) to the data, yielding a best-fit value of a v = 0.48 
ms. This is a factor of ~ 2 less than expected, probably 
owing to the simplicity of the model. Note that the decay 
of 77 observed here is unproblematic for our experiment 
in Ref. [l6| where we took data only for £ s tore < 300 
/is. Furthermore, the decay of a observed in Fig. 2] is 
much slower than the decay of 77 observed in Fig. [6] This 
means that at long storage times, only very little light is 
retrieved but it still has the correct polarization. 

Achieving a much slower decay of 77 during storage 
would be possible when using co-propagating or almost 
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FIG. 6: Decay of the write-read efficiency r\ as a function 
of storage time in a pure BEC. The line shows a Gaussian fit 
according to Eq. ()13p . yielding a best-fit value a v = 0.48 ms 
for the e -1,/2 time of rj. The dominant mechanism that sets 
this time scale is given by spatial filtering caused by the single- 
mode fiber combined with photon recoil incurred during the 
storage process. 



FIG. 7: (Color online) Decay of the normalized write-read ef- 
ficiency r\ in the presence of a noticeable uncondensed fraction 
of the gas. Thermal motion causes a rapid initial decay of r\. 
The decay settles to the long-lived level of rj that is caused by 
the BEC fraction. The data sets were taken for different BEC 
fractions and, correspondingly, for different temperatures. 



co-propagating beams, as mentioned above. This would 
be incompatible with the present atomic level scheme. 
But a conversion of the polarization qubit into different 
wave vectors of the control light, as in Ref. [l2(, could 
solve this problem. However, our experiments in Ref. 
[HI would not immediately profit from a slower decay of 
r\ for two reasons. First, very long storage times would 
drastically slow down the rate at which write-read cycles 
can be repeated, which would result in an unrealistically 
long data acquisition time for the complete experiment, 
due to low count rates. Second, in the setup used in Ref. 
[l6| , the limiting factor when extending the storage time 
was the deterioration of the fidelity due to magnetic field 
noise acting on the single atom inside the high-finesse 
cavity. 

To illustrate how our experiments profit from the use 
of a BEC, we deliberately prepare an atomic gas with 
a noticeable uncondensed fraction. Figure [Jj shows that 
in this case r)(t s toie) decays much faster. More specif- 
ically, rj is the sum of two contributions, one from the 
BEC and one from the uncondensed fraction. These two 
contributions to rj decay on quite different time scales. 
On the time scale shown in Fig. [TJthe contribution of the 
uncondensed fraction decays almost completely, whereas 
the contribution of the BEC is essentially constant. The 
overall decay of r\ is sensitive to the first-order spatial co- 
herence function of the gas [H, HH] . A bimodal decay of 
the first-order coherence similar to Fig. [7] was previously 
observed in Ref. (36j using a different technique. 

Comparing the data in Fig. [7] to the temperatures ex- 
tracted from the size of the uncondensed fraction in time- 
of- flight images, we confirm that AdB/vrci can be used as 
a coarse estimate for the time scale of the decay of rj 
caused by the uncondensed fraction of the gas. Here, 
v Te i = h(k p — k c )/m is the velocity of the F — 2 atoms 



relative to the F — 1 atoms during i s torc, k p is the probe- 
light wave-vector, XdB — y^2Trh 2 /rnksT is the thermal 
de-Broglie wavelength, T is the temperature, and is 
the Boltzmann constant. The data in Fig. [7J were taken 
with a waist of the probe beam of w = 30 /im and after 
removing the filter cavity as well as the single-mode fiber 
between the BEC and the detector. Without these mod- 
ifications, rj would be sensitive only to the central region 
of the gas, where the uncondensed fraction contributes 
less. 



V. CONCLUSION 

To conclude, we characterized and optimized the BEC 
as a quantum memory and showed that a write-read 
efficiency above 50 % can be reached. Its decay over 
storage time results from the differential photon recoil 
in the Raman transfer combined with spatial filtering 
of the retrieved light. This could be mitigated using 
co-propagating light beams. We also showed that the 
mapping between photonic and atomic qubit has an im- 
pressive average process fidelity. Its decay over storage 
time is due to magnetic-field noise and is suppressed by 
appropriate techniques. 
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Appendix A: Modeling the efficiency 

In this appendix, we derive a simple estimate for the 
write-read efficiency r\. This estimate was used for fit- 
ting to the experimental data in Fig. [5J To set the stage 
for this calculation, Sec. lA ll brieflv summarizes the theo- 
retical background from the literature. We then develop 
a simple, largely analytic model for 7/ in Sec. IA 21 In 
Sec. IA 31 we discuss why a large single-photon detuning 
reduces r\. 



1. Theoretical background 

Our notation largely follows Ref. Q, except for the 
sign of all detunings. We denote the probe and control 
Rabi frequencies as fi p and f2 c and the corresponding 
detunings as A p — uj p — ut p . rcs and A c = lu c — w c ,rcs, 
where uj v and lo c are the angular frequencies of the light 
fields and uJ Pties and w c . rcs are the corresponding atomic 
resonances, respectively. We consider the regime of small 
probe intensity and neglect dephasing. The propagation 
of the probe light can then be described by the linear 
susceptibility 2] 



2S 2 T 



X = Xo 



W c - 46 2 (A C + 6 2 ) - 2i6 2 T 



with 



Xo 



"gr"c 

uj p T 



and the group index [2 



ill 



gac. 



(A2) 



(A3) 



Here S 2 = A p — A c is the two-photon detuning, T is the 
total decay rate of the excited state, T p is the partial 
decay rate of the excited state into the ground state in- 
volved in the probe transition, c is the vacuum speed of 
light, g is the particle density, and a — 3A p /2ir is the 
resonant light scattering cross section for a cycling tran- 
sition at wavelength A p = 2irc/io p . 

If the two-photon detuning 6 2 is small, then one can 
expand Re(x) and Im(x) to lowest non-vanishing order 
in S 2 , yielding 



2r 

Mx) = X0^S 2 + O(Sl), 



Im(x) 



Xo 1 ^S 2 



O(Sl). 



(A4a) 
(A4b) 



Note that this is independent of A c . The group velocity 
for probe light can be calculated from Eq. (|A4aP , yielding 



1 + n gI 



(A5) 



We consider a medium that extends from z = to z = L 
with a density g(z) which varies along the propagation 
direction z of the probe beam. If £l c is constant in time, 
the pulse delay after propagation through the complete 
medium follows immediately from Eq. (|A5|i . yielding Q 



dz- 



r 



d p (L), 



where 



f L r 

d p (L) = / dz-^-ag 
Jo 1 



(A6) 



(A7) 



denotes the optical depth seen by the probe light. 

Irreversible absorption of the probe light inside the 
medium can be a serious issue. This is avoided if all 
relevant frequency components of the probe light are 
close to the two-photon resonance, ^2 = 0. Under this 
condition, Eq. (|A4bl) yields a Gaussian EIT intensity 
transmission window in frequency space with rms width 

Ctrans = Awtxans/V^ with [2| 



(L) 



(A8) 



(Al) If Au p denotes the typical width of the frequency spec- 



trum of the probe pulse, then the condition for small 
absorption reads Au p -C Acj t rans(-k)- Due to the Fourier 
limit, the typical duration t p of the pulse is related to 
Auj p by t p Auj p ~ 1. Combining this with Eqs. (IA6|) and 
(|A8|) , the condition for small absorption can be rewritten 
as 



r d (L) 



« \jd P {L). 



(A9) 



Obviously, fully compressing the pulse longitudinally into 
the medium requires r < j(L)/r p > 1. If one simultaneously 
wants to avoid absorption, then according to Eq. (IA9|) 
one needs d p {L) 3> 1. The requirement of large opti- 
cal depth is independent of Q c . But in the experiment, 
il c must be adapted to the values of t p and d p (L), as 
discussed now. 



2. Simple estimate for the efficiency 

While various numerical models for a thorough analy- 
sis of 77 have been published, we find it useful to comple- 
ment these elaborate models with a much simpler model 
that captures only part of the physics but gives a quick 
estimate for the efficiency. Our model assumes that the 
probe pulse is Gaussian in the time domain. Using Eq. 
(|A4[) . the intensity of the propagating probe pulse can be 
approximated as 

I(t, z) = ib(*)exp ("2^) (* - Td ^ - Z ~S) ' ( A1 °) 
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where 1$ is the peak intensity, r p is the temporal rms 
width of the intensity, and r d is given in Eq. (|A6I) . Note 
that Eq. (|A10j) is valid for all z. It simplifies to Eq. ([12]) 
for those values of z where the probe light has not yet 
entered the medium. 

Our simple model separately addresses the issues of in- 
sufficient pulse compression and irreversible absorption. 
First, we ignore irreversible absorption. This makes Iq 
and t p independent of z. We assume that Q c is switched 
off at a time to and assume that the fraction of the light 
that is inside the medium at this moment is stored and 
subsequently retrieved. This yields a write-read efficiency 



1 



/comp 



erf 



to 



V2t 



-erf 



t - T d (L) - L/c 
V2r p 



(All) 



with the error function erf(x) = (2/y / 7r) J Q X duexp(—u 2 ). 
This result for the efficiency quantifies how well the pulse 
is compressed into the medium. Note that n gr 3> 1 im- 
plies T d (L) > L/c. 

Second, we turn to irreversible absorption. Here, we 
consider a situation in which Q c is constant in time, im- 
plying that no storage takes place. The fraction Straus of 
the pulse energy that is transmitted through the medium 
is calculated easily in the frequency domain. Based on 
Eqs. (IA4bp and (|A10I) . we obtain 



Vtta 



1 



(TpAw trans (L)) 5 



-1/2 



(A12) 



with Awtrans from Eq. (|A8|) . This result for the efficiency 
expresses the issue of irreversible absorption. 

To obtain a simple estimate for the overall efficiency, 
which must take both effects into account, we simply 
multiply the two efficiencies from Eqs. (IA11[) and (|A12I) . 
For the atomic probe transition used in our experiment, 
we have T p = T/12. Combination with the atom number 
and trap frequencies quoted in Sec. Ill El yields a peak 
value of d p (L) = 127 at x = y = 0. Results for this 
optical depth and typical parameters of our experimental 
pulses are shown in Fig. [5J 

In the model, we explored rj as a function of the two- 
dimensional parameter space spanned by Q c and to- Fig. 
[8J shows the dependence on Q c only for that value of 
to, for which the global maximum of r\ is reached. A 
modification of the pulse duration t p would require a re- 
optimization of Q c and to- Considering Eqs. (|A6[) and 
(fA"8j) . one finds that Eqs. (|A11|) and (|A12I) remain un- 
changed if the scalings Q c oc l/Vra and to oc t p are used. 
As a consequence, the maximum value of ri(fl c ,to) is in- 
sensitive to a change in r d . 

The transverse inhomogeneity of the BEC can be ac- 
counted for by calculating a weighted average of rj, with 
the transverse profile of the probe light intensity as a 
weighting function 



^comp 
^trans 

Icomp^trans 
average 




10 15 20 25 30 35 40 



Q c /27t (MHz) 



dxdy- 



-2(x 2 +y 2 )/w'- 



r}(d p (x,y)). (A13) 



FIG. 8: (Color online) Simple theoretical estimate for the 
write-read efficiency rj vs. Rabi frequency of the control light 
f2 c . The dashed, dotted, and solid lines show 77 comp , 7/trans, 
and ?7com P X ?7trans, respectively. Parameters are t p — 94 ns, 
to = 230 ns, 1/r = 26 ns, and d p (L) — 127. The pronounced 
maximum in the solid line arises because for small Q c , the 
EIT transmission window in frequency space is too narrow 
for the incoming probe pulse, whereas for large fi c , the group 
velocity is not sufficiently reduced to spatially compress the 
complete pulse into the BEC. The dash-dotted line shows the 
result of averaging over the transverse profile of the BEC. 



We assume a Thomas-Fermi parabola for g(x, y, z) with 
Thomas- Fermi radii R x , R y , and R z . Calculation of 
d p (x, y) by analytic integration over z is straightforward. 
After a transformation to new coordinates (p, ip) with 
x = R x p cos ip and y — R y psinip, the integral over ip can 
also be solved analytically. The remaining integral over 
p is easily computed numerically. The result is shown as 
a dash-dotted line in Fig. [5J This line predicts a maxi- 
mum of r\ ~ 60 % at f2 c = 27r x 15 MHz. For this value 
of f2 c and for x = y — 0, the theory yields estimated 
values of xo — 0.5, n gr = 5 x 10 6 , Td(L) = 550 ns, and 

Awtrans(i) = 2?T X 3.3 MHz. 

The absorption represents a filter in frequency space. 
Due to the Fourier limit, this causes an increase of the 
temporal pulse width r p for increasing z, thus complicat- 
ing a more rigorous calculation of rfcomp ■ We can overes- 
timate the effect of this increase of t p (z) when using the 
constant value t p (L) instead of t p (0) in calculating rycomp- 
We find that for the parameters of Fig. [8j this has little 
effect. Our model also neglects that t p (z) should increase 
due to dispersion caused by Re(x). This increase is given 
by if{z) = 7-2(0) + [T d (z)/T p (0)u} p ] 2 , which is negligible 
for the parameters of our experiment. 

The simple model developed here neglects that the spa- 
tial cutting during storage broadens the frequency spec- 
trum of the pulse, thus increasing the absorption after 
retrieval. The model also neglects that the transverse in- 
homogeneity of the medium might cause a deformation 
of the wavefronts, resulting in effects such as focussing of 
the probe beam. In addition, this model is fully based 
on Eq. (|A4p . instead of Eq. (|A1[) . If the probe pulse is so 
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FIG. 9: (Color online) Linear susceptibility x as a function 
of the two-photon detuning 82- The solid and dotted lines 
show the predictions of Eq. (IA1|) for the real and imaginary 
parts of x at n c = 3.31" = 2tt x 20 MHz. (a) At the single 
photon resonance, A c = 0. (b) At a single photon detuning of 
A c = 11.4T = 27r x 70 MHz. Dash-dotted lines (red) show the 
lowest-order approximations of Eq. I|A4(I . The range where 
these approximations become poor sets an upper bound for 
the frequency range which is useful for light storage. This 
range is much narrower in (b) than in (a). This explains the 
reduction of the efficiency observed at A c = 2n x 70 MHz. 
Note that the scales on the horizontal axes differ by one order 
of magnitude. 

broad in frequency space that Eq. (|A4[) is not a good ap- 
proximation, then the applicability of the model is ques- 
tionable. Nonetheless, the prediction of the model agrees 
reasonably well with our experimental data in Fig. [SJ 

3. Reduced efficiency at large single-photon 
detuning 

To understand the physical origin of the reduction of 
r\ at large single-photon detuning, we investigate the fre- 



quency range that is useful for storing light. An up- 
per limit for this frequency range is set by the frequency 
range within which Eq. (|A4|) is a good approximation. 
Using Eq. (|A1|) , one can easily show that the maxima of 
Im(x) are located at 



h= 2 (-A c ±VAf+fif). (A14) 



The maximum nearest to 82 = clearly sets an overop- 
timistic upper bound for the useful frequency range for 
light storage. For |A C | ^> fi c , the nearest maximum lies 
at S2 ~ Og/4A c . Comparison with Eq. (IA8|) shows that 
the useful frequency range is much narrower than Awtrans 
unless 4|A c |/r^/d^ <C 1. Our experiment is operated at 
4A c /r-\/d^ 4 so that this issue is obviously relevant. 
This reduction of the useful frequency range is to be con- 
trasted with the delay Td in Eq. (|A6[) which is indepen- 
dent of A c . As a result, the overall efficiency is reduced. 

To further illustrate this point, we show the depen- 
dence of x on ^2 in Fig. [9] Parts (a) and (b) correspond 
to A c = and A c = 2n x 70 MHz, respectively. Note the 
different scales on the horizontal axes. This figure clearly 
illustrates that the frequency range over which Eq. (IA4|) 
is a good approximation differs drastically between the 
two cases. 

The light pulses that we store at A c = 2n x 70 MHz 
have the same spectral width as for A c = 0. For A c = 
2ir x 70 MHz a considerable part of the frequency compo- 
nents of the light therefore samples the frequency range 
where Eq. (|A4[) is not a good approximation. Fig. [9fb) 
shows that for components with negative 62, the value 
of dRe(x) I dS2 is reduced, resulting in a faster group ve- 
locity, which is disadvantageous. On the other hand, for 
positive 62 absorption can be substantial and dRe(x) / (IS2 
can even change sign, thus not creating slow light. These 
problems qualitatively explain the reduced write-read ef- 
ficiency that we observe experimentally for A c — 2ir x 70 
MHz. 
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